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Quasi 1-D analysis of a circular, compressible, turbulent jet laden 

with water droplets 


September 27, 2001 


Abstract 

Recent experimental studies indicate that presence of small amount of liquid droplets reduces 
the Overall Sound Pressure Level (OASPL) of a jet. Present study is aimed at numerically in- 
vestigating the effect of liquid particles on the overall flow quantities of a heated, compressible 
round jet. The jet is assumed perfectly expanded. A quasi- ID model was developed for this pur- 
pose which uses area- averaged quantities that satisfy integral conservation equations. Special 
attention is given to represent the early development region since it is acoustically important. 
Approximate velocity and temperature profiles were assumed in this region to evaluate entrain- 
ment rate. Experimental correlations were used to obtain spreading rate of shear layer. The 
base flow thus obtained is then laden with water droplets at the exit of the nozzle. Mass, mo- 
mentum and energy coupling between the two phases is represented using empirical relations. 
Droplet size and mass loading are varied to observe their effect on flow variables. 


Nomenclature 


p g \ Bulk mean density of gas phase [kg/m 3 J 
u g : Bulk mean velocity of gas phase [m/s] 
m p : Mass of one particle [kg] 

m e : Entrained mass per unit length in streamwise direction [kg/m-s] 
P: Pressure [Pa] 

T g : Bulk mean temperature of the gas phase [K] 

Tp : Temperature of liquid particles [K] 

D: Nozzle diameter [m] 

Up : Particle velocity [m/s] 

n: Particle flux through any cross section of the jet [s~*] 

A: Effective area of the jet [m 2 ] 

Cpg : Specific heat at constant pressure of the gaseous phase [J/kg-K] 
L: Latent heat of vaporization of water [J/kg] 

Cp: Specific heat of liquid water [J/kg-K] 

R: Gas constant per unit mass [J/kg-K] 
h: Convective heat transfer coefficient [J/m 2 sK] 
ho' Mass diffusivity [kg/m 2 s] 

D m : Diffusion constant [m 2 / s] 

C2: Non-dimensional entrainment rate 
Cp: Drag coefficient 
a: Entrainment coefficient 
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M c '. Convective Mach number 
a: Radius of droplet [m] ; 

Speed of sound [m/s] 

5gg: 99% velocity thickness [m] 

81 / 2 : 1/2 velocity thickness [m] 

8 „: Vorticity thickness [m] 

8 V is : Visual thickness [m] 

Spit * Pitot tube thickness [m] 

M: Momentum flux [kg — m/s 2 ] 
m: Mass flux [kg/s] 

r: Radial distance from jet centerline [m] 

Re: Reynolds number 

Pr: Prandtl number 

Sc: Schmidt number 

Bi: Biot number 

Sh: Sherwood number 

Nu: Nusselt number 


Subscripts 

g: Gaseous phase 
p: Particle 
e: Entrainment 
w: Liquid water 
v: Water vapor 
a, air: Air 
sat: Saturated 
0: Jet exit 
c: Centerline 
00 : Ambient 


1 Introduction 

Supersonic jet noise can be classified into three distinct components- shock associated noise, 
screech tones, and turbulent mixing noise (Krothapalli and Washington, 1998). The turbulent 
mixing noise is the most dominant contributor to the OASPL. A novel approach to jet noise 
reduction is to use water injection to minimize the turbulent mixing noise. The goal here is to 
use minimal mass loading of the gas phase with water (less than 10%) to achieve substantial 
reduction in OASPL. Higher mass loading is impractical as the aircraft has to carry the water 
to be injected into the jet. 

A 1-D model of the jet was developed. It uses continuity, momentum and energy equations to 
march variables in the streamwise direction. Area^ integrated form of the conservation equations 
are used. Variables considered are area-averaged flow variables. The assumption that these 
equations represent the actual system to a fairly good accuracy can be justified by the following 
analysis. 
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1.1 Integral form of RANS equations 


Since algebra can grow quite tedious with cylindrical coordinates, we will consider a plane jet 
for illustration purpose. Also, to simplify the analysis, we assume incompressible flow. First, 
consider continuity equation for turbulent flow. 


Integrating over y, 


but v\h — and 

hence, 


du dv 
dx dy 


— 0 


£(/ ud y) +v tf = ° 

£{J 0 


v|o = 0 by symmetry. 


A 

dx 



rn e 

P 


Now consider the x-momentum equation in conservative form- 


a) 


du 2 duv d /—*\ d (— \ 1 dP 

Integrating over area perpendicular to the x-axis, 


fd 2 u d 2 u \ 
\9x 2 + dy 2 ) 


W« ^YY)l + Wo t *) + ( uV )| ud >) + ^ 

Last term on each side, and the second term on the LHS are zero, when evaluated at y=0, 

H. Since mean pressure is assumed to be constant everywhere, dP/dx = 0. Also, for a turbulent 
jet, viscous stresses are unimportant. So, simplified form of the above equation becomes: 


H 


A 

dx 




We know that v! is less than roughly 10% of u in magnitude. Hence second term in the 
bracket is even smaller than the first term. The above equation is then approximately equivalent 
to 


Wo ^Y 0 <2> 

This is the integral form of the momentum equation. Energy equation can also be treated 
similarly. We will apply these equations to the present case of compressible jet using mean 
variables. In order to completely specify the problem, one needs to specify the entrainment 
mass rate distribution. To obtain this, experimental data about spreading rates of jet and its 
shear layer are used. Such data is presented by Papamoschou and Murakami (2000), Lele and 
Freund (2001) and Dimotakis et al (2000). Experiments related to prediction of entrainment rate 
were done by Hill (1972), and Ricou and Spalding (1961). Relation between different measures 
of spreading rates, namely, vorticity thichness, visual thickness and pitot tube thichness was 
obtained from Papamoschou and Roshko (1988). Experimental and analytical study of jet 
entrainment is also found due to Zaman (1998). Further, DNS study of particle-laden mixing 
layer by Miller and Bellan (1999) was found to be a useful reference for phase coupling relations. 
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2 Problem formulation and governing equations 

A circular jet issuing into quiescent atmosphere is considered. The geometry and jet exit condi- 
tions were those used in the experiments by Krothapalli and Washington (1998). Jet diameter 
is 29 mm. The jet is assumed to be perfectly expanded i.e., pressure everywhere outside the 
nozzle is equal to atmospheric pressure and the nozzle is operating under design conditions. The 
jet Mach number is varied as a parameter. Mach number at the exit determines the exit and 
stagnation conditions through the requirements of isentropy. Jet Mach number was kept 1.44 
in this study. Stagnation temperature was set to 600 K. Liquid water was introduced in the 
form of tiny spherical droplets at the jet exit at ambient temperature. Ambient temperature 
was taken as 298 K. Initial velocity of the droplets at exit was 50 m/s. ”n” number of water 
droplets of instantaneous mass m p cross any given cross section of the jet per unit time. We can 
extend our analysis to accomodate different droplet sizes and corresponding number densities. 

As discussed earlier, approximate integral equations will be used for marching. These 1-D 
equations with the streamwise coordinate x as the independent variable are written in control 
volume form. Entrainment rate is calculated simultaneously (discussed in detail later), and the 
equations are spatially marched. But before doing that, we define the equivalent plug profiles 
for velocity, density and temperature and equivalent area of the jet in such a way that these 
quantities confirm to the conservation equations and constitutive relation. These variables are 
then used as dependent variables which are solved for. Definition of such variables is illustrated 
in figure (1) below. 



Figure 1: Defining area averaged variables 

Continuity equation: 


Tvm p 

+ (PgAUg) 

+m e 

Ax = ( pgAu g ) 

I 7lT7l p 


X 

X 

X 

x+Ax 


dx (p 9 Au 9 ) + &E ( nm p ) — m e (3) 

Momentum balance: 

Momentum in + Force in the -fve x direction = Momentum out 


{P9 U 1 A )\ 


+ (nm p u p )\ 


= {p g u]A) 


+(nm p u p ) 


x+Ax 


Ix+Ax 
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Net pressure force is zero because pressure is constant everywhere (atmospheric) and the 
control volume is closed. Above equation then becomes, in differential form, 


Energy balance: 



( 4 ) 


( Thermal + KineticEnergy) 


= ( Thermal + KineticEnergy) 


tin 


out 


(PgUgACpgTg ) 


+(nm p C p T p ) 


+{±PgU 3 g A) 


— (pgllgACpgTg) 


x+Ax 


+(nmpCpT p ) 


x+Ax 


+( nm Plf) 

+ (lp g U 3 g A) 


X ~h ni^Cyp^air'Eamb^X 
+{nmp^f) | 


x+Ax 


Ix+Ax 


^ x (p g u gAC P gT g ) + — (nm p C p T p ) + = m e C p <air T aTnb (5) 

In addition to these conservation equations, we also have relations describing the coupling 
between the two phases corresponding to mass, momentum and heat transfer. Rate of mass 
transfer between a water droplet and the surrounding gas is determined by the difference between 
concentrations of water vapor at the surface and away from the surface. It is governed by the 
relation: 


= -{Area)h D ln{\ + B m ) 

^Ef-Up = -4Tra?h D ln(l + B m ) 

^ = -=^-h D ln(l + B m ) ( 6 ) 

where B m = » X&at and X v being mass fractions of vapor under saturated and actual 

conditions respectively. /i£> is mass diffusivity of water vapor in air in kg/m 2 s . 

X 3a t and X v are calculated as follows: 

If m v is the mass flux of water vapor through any cross section of the jet, then 

_ m v At m v 
AUgAt Au g 

Since water vapor is generated only due to evaporation of droplets, 

TTiy — n(m po Trip ) 

Water vapor density under saturated conditions is 

_ Pi sat 
PS<>t ~ RvTg 

with Ry being the gas constant for water vapor and P sat the saturation pressure given by 
(Miller and Bell an, 1999) 

Psat = Pat m exp{jj-(P-X)} 
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Here Tb is the boiling temperature of water at atmospheric pressure, L is the latent heat of 
vaporization of water (assumed constant), and T p is particle temperature. Now, patial pressure 
of dry air under saturated condition is P a ,sat = Patm — Psat • So, the density of dry air in 
saturated air becomes 


Passat 


P z,sat 
RairPg 


Then, the total density of gaseous phase will be 


Ptot,sat — Passat “1" Psat 


The mass fractions now become 


y Psat 

^ sat = 

Ptot f sat 


CLTld Xy 


Pv 

Pgas 


Mass diffusivity is given by 


hj^ — Sh pg 


Pm 

(2o) 


Where 

Sh = Sherwood number 

An^Diffusivity of water vapor in air = 2.93 x 10 -5 m 2 /s 
2a=Length scale (droplet diameter) 

Sherwood number is dependent on Reynolds number based on the empirical relation 


Sh = 2 + 0.552 Sc 1/3 Re 1/2 (8) 

Here, Sc is Schmidt number and its value is 0.608 for water vapor in air. Reynolds number 
is defined based on the particle diameter and slip velocity. 

Re= Pj l"j,-“pl(2°) (9| 

Pg 

Momentum coupling relation is based on experimentally measured values of drag coefficient, 
Cd, for a range of Reynolds number. Drag force is 

F d = C D ~Pg(u g - u p ) 2 (na 2 ) (10) 

But Fd — mpdu p /dt = m p Updu p /dx. Hence, 

dUp 1 1 . \ 9 / 2 \ \ 

~dx = m p Up^ D 2^ 9 ^ U9 ~ ^ Ka ) s ^ n ^ u 9 ~ u p ) (11) 

Variation of Cd with Re used here is shown below (Donley, 1991). 

Heat transfer between the two phases is governed by the temperature difference between the 
particle surface and the surrounding fluid through the relation 

h(T p - r 9 )W = ^[L + C p (T evap - T p )) - m p C p ^ 

Substituting At = A x/u p , 


h(T P - T g ) 


47ra 2 

Up 


^[L + C p (T evap — T p )] - mpCp ^ 


( 12 ) 
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Figure 2: Variation of drag coefficient with Reynolds number 


Where L is the latent heat of vaporization of water. Its value is 2.5 x 10 6 J /kg. Here it is 
assumed that the temperature inside a droplet is approximately constant i.e. temperature gra- 
dients inside a particle are small. This is justified when particle size is small and its conductivity 
is high. A plot of particle Biot number evolution (Fig. (3)) for the particle size of 26 microns 
shows that Bi = hd/6h w <0.1 and above assumption is valid. 



Figure 3: Evolution of Biot number 
To compute convective heat transfer coefficient we use the following relation 

Nu = 2 + 0.552Pr 1/3 Re 1/2 (13) 

Now, changes in thermodynamic properties of the gaseous mixture are functions of its com- 
position i.e. fraction of water vapor present. Variations in these properties are governed by the 
requirement that, in an adiabatic mixing process, total energy is conserved. We define the bulk 
properties of the gas phase in such a way that these properties, when used to obtain energy, 
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give the same total energy as the sum of the energies of the individual gaseous components, 
assuming a hypothetical adiabatic mixing process. 

Applying the above statement to the control volume, 


(pgUgACpg) 

T 

— ( PgUgACpg ) 

T 

+77i e 

AxC P: airT 


x-f-Ax 

x+Ax 

X 

x+Ax 

X 


+n- 


dm T 


x+Ax 


dx 




x+Ax 


'^j—\P9 u gAC’pg) = '™ , eC'p,air H" ^ Cp,v 

Same equation can also be obtained by mass- averaging the component Cp values. 
It should be noted that partial pressures of gases add to produce total pressure. 


(14) 


P = P a + P v 


Also, densities add. 


P — {pa.Ra “ 1 “ PvRv)T g 


(15) 


Pg — Pa “ 1 “ Pv ( 16 ) 

Thus, we have a set of the above equations and the following unknowns: u g ,p g , T g , C pg , u v , m p , T p , A. 


3 Solution procedure 

Solution was marched streamwise using Euler explicit method. This required derivatives of some of the 
dependent variables at every step. It was observed that this produces an ill-conditioned system of equations 
that tends to diverge. As a remedy, this system was diagonalized by defining two more variables related to 
the gaseous phase, namely, mass flow rate and enthalpy of the gas. 


Trig pgtLgA 

hp = CpgTg 

Henceforth, we adopt a notation in which a dot (.) above a variable indicates derivative with respect to 
x. Variables pertaining to particles, u p , T p , m p can be easily decoupled using the coupling relations in the 
following way: 

Equation (6), (11) ,(12) respectively give: 


rn P = -^h D ln( 1+B m ) 

(17) 

u 'p = ^jz; c D%Pg(ug - u p ) 2 na 2 sign{u g - u p ) 

(18) 

Tp — nipCp Tg) Up ~ ripiL + Cp{T eva p “ T p )]| 

(19) 

Variables corresponding to the gas phase are then decoupled as follows: 
Equation (3) gives: 


rfig — Trig TiTYip 

(20) 
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Equation (4) is: 


£(m g Ug) + £{nm p u p ) = 0 

TTlgUg -f Ugrilg + TlfapUp -f UpTYlp) = 0 
TTlgUg UgTTlg Tl^TTlpUp “f~ lipTflp} 

Using eq.(20), we get 


TTlgUg = —u g m e -f nUgThp — nrfipUp — 7iu p 7n p 
7n g u g — —u g m e + n(u g — Up)m p — nu p m p 
u g = ^ [—u g m e + n(u g - u p )m p - nu p 7n p } 


Energy balance equation was: 


A 

dx 


(m.gh g ) + 



— TTl e h e 


AVhere /l e — ^p^air^amb 
Rearranging, 


h g {^e^e TTlgUg TTlgUgUg 71 [?7ip ( Cp Tp ~f~ “h TTlpCpTp -J - T/ZpUptipj } 

771 g & 2 

Equation (14) was: 

’d^iPg^gAC'pg) = 7n e Cp^air ”1“ 71 

This gives: 


(21) 


(22) 


m gCp,g — 7Tl e Cp i air “f" TlTTlpCp ^ ^ , g^'p,g (23) 

Equations (17) to (23) form a diagonal system of linear algebraic equations in u p , T p , m p , u g , m g , h gi C pg . 
After advancing every time step, all other variables are updated as: 


Tg — hg/Cp,g 

Pg — Pa ^ pv 


Where p a =^ and Pv = ^ 


(24) 


^aif — { Jq TTtgdoC 

TTlvapor ~ 7l{rTlp 0 /Up) (25) 

This system of equations can be regarded as an implicit system in area A. We first linearly extrapolate 
the area to the next node, evaluate all other variables, and then correct the area using the following set of 
equations: 
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Pg Pa “I" Pxu 


A = 


PgUg 


In this way, we can iterate to get the converged value of the area. Here it was observed from the following 
plot of order of magnitude of relative arror in area as a function of number of iterations (Fig. (4)) that only 
two iterations suffice to achieve desired accuracy. Using any further iterations would be an improper use of 
computational resources. 



Figure 4: Variation of relative error in area with number of sub-iterations 


4 Calculation of entrainment rate 

In the above solution procedure it was assumed that the entrainment rate was supplied at each marching 
step. This section focuses on obtaining the rate of mass entrainment. Developing region and self-similar 
region are considered separately. 


4.1 Fully developed region 

Consider incompressible flow. In the self-similar region, the jet spreading rate in terms of the half- velocity 
thickness is given by the expression (Lele and Freund) 


^ « 0.095-^ 
D D 


„ 2 \ -2 


The self-similar velocity profile is 

f=K)' 

where r; = r/5 sc and S sc is the scaling thickness equal to 0.38845 j/ 2 (White). Momemtum flux is 


(26) 


(27) 
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M = pU?5lf~2nT,£dr, 

= pU?5lf ~2Tt V (l + ^ 4 dr, 
= P U?5l( 4.1888) 


Hence, 


(28) 


n = J_ / M 

c 5 SC Y 4.1888p 


(29) 


An entrainment hypothesis was made which states that the turbulent fluctuating velocity that causes 
entrainment is proportional to the jet centerline velocity. 


77i e — aU c p s 27rS sc (30) 

where p 3 is the ambient density and a is the entrainment constant to be determined using the experimental 
data related to the jet spreading rate. This gives, 


^ = aU c p2n5 sc 

=> jt^seP/c) 00 27T7?^1 + <fyj = ap s U c 2n8 sc 

■* fc 2 ’ 

But 5 SC = 0.38845X/2 = 0.0369x. This gives, 


a = 0.0738 

We now define the non-dimensional entrainment rate as C 2 by 



m Q is the mass flow rate at nozzle exit and p 0 is fluid density at nozzle exit. Here, p Q = p 8 — p since we 
have assumed incompressible flow. 

° 2 = pU c nD 2 /4 (pi ) 2 ™\j 4.1888p P (32) 

But M = jD 2 pU%. The above expression then becomes 


C 2 


^C 2 = 0.256 


pU c nD 2 /A 


OO 


1/2 


2 -kolDUc 


j 7 r 

y 4 4. 


1888p 


= 4a 


4.1888 


(33) 
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Ricou and Spalding (1961) state that C 2 ranges from 0.248 to 0.456. Experimental study due to Hill (1972) 
suggests the value 0.32 for C 2 . Zaman (1998) predicts a similar value, 0.316. Since C 2 varies considerably 
in different studies, we adopt a representative value of 0.32. 

4.2 Early development region - Compressibility effects 

The developing region requires special attention as it is important from acoustics point of view. Compress- 
ibility effects in this region may be significant. They are incorporated into the present analysis. As the 
jet issues out of the nozzle, velocity profile is a plug profile. As the vortex sheet at the edge of the profile 
diffuses, the annular shear layer spreads. Spreading rate depends on the convective Mach number defined as 

U c 

M c = — — ± (34) 

0>e t (lamb 

Where U c = Jet exit velocity 
a e — Speed of sound inside the jet at exit 
c iamb = Speed of sound in atmosphere 

Papamoschou and Murakami (2000) give an expression for spreading rate of pitot tube thickness: 

S^t = 0.14(1 + v/s)(0.23 + 0.77exp(-3.5M c 2 )] (35) 

where S — pamb/ Pexit 

In the early region, sinusoidal velocity profiles were assumed in the shear layer. These profiles are: 



C p T — — — + Au + B (37) 

Applying the boundary conditions: 

At r = 0, u = J7 C , T = T c 
At r = 5, u = 0, T = T c 
we get 



(38) 

(39) 
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Using these profles and the spreading rate data, momentum conservation can be applied to obtain thick- 
ness of the potential core. This completely describes velocity and temperature profiles. This information 
can then be used to compute entrainment rate using mass conservation. If M is the momentum flux at the 
exit of the jet then at every cross section along the jet, momentum flux remains M. This gives 


rVc ry c +6 

I pu^nrdr + / pu 2 2nrdr = M 
JO Jy c 


(40) 


Denoting (r — y c — by (5 and using the ideal gas law, we obtain: 


pcU^yl + fj f -gfu 2 + y c + %d /3 

y 2 c + $i&{ [A tv] v4 + A t (f + i) }= ¥ 


(41) 


Here we have used the fact that 


D 2 

M = 7r——p c U i 


This is a quadratic in y c (x) and can be solved. Integrals appearing above are evaluated numerically and 
Spit(x) is supplied. 8 is computed from the pitot tube thickness by using the equations (Papamoschou and 
Roshko, 1988): 


X . ~ °Pit 

o m3 ~ 0 8 

8(jj ~ 0 . 5^3 


(42) 

(43) 


Having obtained velocity and temperature profiles completely, continuity equation can be used to compute 
entrainment rate. We have, 


m e = 


dm d f y * +s „ , 

= — / puznrar 
dx Jq 


dx 


Using Leibnitz rule, 


= Jo 


Trip 


dr + ( pu2nr ) 


l/ c +<5 


■^(Vc + 5) - (pu2nr) 




Last two terms on the RHS are zero. 

r Vc +6 Q 


m e 


[ Vc+ d d(3 f Vc d 

=1. &**.**+}. dff lpu2 ’ T) 


. dx 


Second integral on the RHS is zero because both density and velocity are constant inside the potential 
core. 

Changing variable of integration to /?, 


Now, 


771 e 



r 2 'w i,u) ^ dl3 


(43) 
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( 44 ) 


“ P d(5 + U dfi 

— i u d£dTdu 

Pd(3 ^ u dT du dp 


But P= jfr giving us: 


Also, T = [— \ + Au + B\/C p . Hence, 


dp _ P 
dT ~ RT 2 


Where 


Also, 


% = (A-u)/C p 


Hri = 


d0 


p + U 


( 

l RT 2 J 


L±Z*1 

C p 


du 

dp 


du _ d 
7p = df3 


y(l - sinfi) 


U c 

= ~cos/? 


it = _ F r,S ' + fiVc# - f y'c (46) 

Where dash denotes differentiation with respect to x. Now, m e can be evaluated by numerically com- 
puting the integral in equation (43). As the shear layer grows further, it starts ” feeling” the effect of its 
circular shape and the growth rate no longer remains constant. In the present study, the above analysis was 
applied upto a distance of twice the jet diameter. Beyond 13D, the constant value of C 2 corresponding to the 
fully developed flow is used. The middle region is interpolated using a cubic spline in the non-dimensional 
coordinates C2 and x/D. Resultant variation of C2 is plotted in figure ( 5 ). 


5 Results and discussion 

To generate data and highlight the effect of liquid droplets on the jet flow, we keep the geometry and exit 
condition of the jet fixed. Jet diameter is 29mm, jet exit Mach number is set to 1.44, stagnation temperature 
is 600 K. This gives jet exit velocity of 594.44 m/s and temperature at the jet exit is 424.1 K. 

First we shall demonstrate certain general trends using the particle size of 26.7 microns. Figure (6) shows 
that the momentum coupling between the two phases causes the particle velocity to tend to follow the gas 
velocity. Heat coupling has a similar effect on the droplet temperature. 

We set the total mass flow rate of water droplets at 10% of the jet exit mass flow rate. We then vary 
the droplet size to observe the droplet behaviour. From figure (7) we see that smaller particles affect the 
velocity history more than the larger particles do. 

Some other important observations are: 
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Figure 6: (a) Velocity and (b) Temperature response of a particle 

• Computations were performed to evaluate the importance of surface energy change during 
atomization. It was found that, if every particle breaks into two, then the increase in 
surface energy is six orders of magnitude less than the kinetic energy of the gaseous phase. 
This suggests that, the mechanism of conversion of turbulent energy into surface energy 
is unimportant here. 

• A simple modification of the formulation was done in order to accomodate different par- 
ticles sizes with the corresponding number densities. 

• Efficacy of Crocco-Busemann relation in conserving total energy was tested. Ass um ing 
velocity profiles and temperature profiles are same, the relative error in total energy was 
found to be 3% or less. After using Crocco-Busemann relation to find temperature profiles, 
this upper bound on the relative error in energy was reduced to 0.5%. 

• Dimensional analysis suggested that the non-dimensional parameter deciding the impor- 
tance of particles in mass and momentum variation was nothing but the ratio of the total 
mass flow of liquid droplets to the mass flow rate of gaseous phase. Also, since temperature 
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xJD 


Figure 7: Effect of droplet size on main flow 


variations for both the phases are approximately of the same order, the same parameter 
also governs energy variation. When this parameter is small, presence of the droplets may 
be neglected completely. 
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